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TECHNICAL NOTE 5157 

METHOD FOR CALCULATION OF COMPRESSIBLE LAMINAR 
BOUNDARY LAYER WITH AXIAL PRESSURE 
GRADIENT AND HEAT TRANSFER 
By Paul A. Libby and Morris Mbrduchow 


SUMMARY 


A rapid and sufficiently accurate method, for most practical purposes, 
of determining laminar -boundary-layer characteristics in flow with a given 
free-stream Mach number and given velocity distribution at the edge of the 
boundary layer is presented. The method can be easily applied to flow with 
zero pressure gradient for any (constant) Prandtl number of the order of 
unity and any given temperature distribution along the wall. Numerical 
examples are given to illustrate the method and the satisfactory accuracy 
obtained. For flow in an axial pressure gradient, the method can be 
applied for a Prandtl number of unity and any given uniform wall temper- 
ature. The methods developed here are based on an application of the 
Karmen integral method to both the momentum and energy equations, in 
conjunction with a sixth-degree velocity profile and a seventh-degree 
stagnation-enthalpy profile. A single boundary-layer thickness and one 
of the coefficients in the thermal profile are the parameters in this 
two -parameter method. 


INTRODUCTION 


The aim of this report Is to present a relatively simple method, 
sufficiently accurate for most practical purposes, of calculating the 
laminar -boundary -layer characteristics in the compressible flow over a 
given object with heat transfer at the wall. 

The method is based on the extension of the Karman-Pohlhausen method 
to sixth-degree velocity profiles and seventh-degree stagnation-enthalpy 
profiles. The use of sixth-degree velocity profiles is in accordance with 
the conclusions of reference 1, wherein it was found that such profiles 
can usually be expected to lead to results of adequate accuracy without 
much increase in computational work. Such profiles have been applied 
with satisfactory results for compressible flow over a flat plate with 
heat transfer (ref. 2) and for compressible flow in an axial pressure 
gradient without heat transfer (refs . 3 and 4 ) . 
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The present analysis consists of two main parts. In the first a 
simple approximate general solution is developed for flow without axial 
pressure gradient, with any given Mach number_ (neglecting hypersonic 
effects) and given distribution of temperature along the wall. The 
Prandtl number is also left arbitrary, although it is req uir ed to be of 
the order of magnitude of unity. The coefficients of specific heat, as 
well as the Prandtl number, are assumed as constant, while the viscosity 
coefficient is assumed to vary linearly with the temperature, the propor- 
tionality factor being chosen in accordance with the Sutherland relation, 
as in references 1, 2, 3, 5, and 6. In the second part of this analysis, 
a simple approximate general solution is developed for flow in a pressure 
gradient with any given constant wall temperature, free-stream Mach num- 
ber, and velocity distribution outside the boundary layer. The Prandtl 
number is, however, assumed to be unity. 

A brief indication of the pertinent literature may be worth while. 

The flow without axial pressure gradient (e.g., flow over a flat plate) 
has, of course, already been treated to a considerable extent. Karmen 
and Tsien (ref. 7) derived solutions for a Prandtl number of 1, constant 
wall temperature, and viscosity coefficient proportional to a power of 
the temperature. Chapman and Rube sin (ref. 3) have obtained exact (numer- 
ical) solutions for variable wall temperature and a linear viscosity- 
temperature relation. Van Driest (ref. 8) has accurately calculated a 
series of curves showing skin friction, heat-transfer coefficients, and 
velocity and temperature profiles for a constant wall temperature, using 
the Sutherland viscosity -temperature relation. Klunker and McLean 
(ref. 9) have calculated the boundary- layer characteristics using the 
actual variation of Prandtl number, specific heat, and viscosity coeffi- 
cient with temperature, while Young and Janssen (ref. 10) and Moore 
(ref. 11) have done likewise with the aid of a differential analyzer, 

Moore taking dissociation into account. 

In each of the above references, the calculations, though exact, are 
fairly tedious, while, except for reference 5, the wall temperature Is 
assumed as constant. The advantage of the method presented in the present 
report is that the calculations can be performed with relatively little 
difficulty for a prescribed wall -temperature distribution, Mach number, 
and Prandtl number. Although the results thus obtained will not be exact, 
they will usually be sufficiently accurate in practice. 

In contrast with the literature on flow over a flat plate, the liter- 
ature on compressible boundary -layer flow in a pressure gradient with heat 
transfer still appears to be limited. A good survey of work done on 
incompressible flow can be found in reference 12. Kalikhman (ref. 13) 
has treated the compressible case by applying the Karman-Pohlhausen method 
to both the momentum and energy equations. Ferrari (ref. l4) has recently 
presented an accurate, but elaborate, method of calculating the boundary- 
layer characteristics for arbitrary Prandtl number and uniform wall 
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temperature. The present analysis has, once again, the advantage of 
presenting a general approximate solution which involves relatively little 
computational difficulty. However, the Prandtl number is here restricted 
to unity. 

The method of analysis used in this investigation merits some comment. 
As in references 13 and 15, stagnation-enthalpy, rather than temperature, 
profiles are used. The use of stagnation-enthalpy profiles leads to 
significant mathematical simplifications, while it can be shown, more- 
over, that for a Prandtl number of 1 the equations used here automatically 
lead to what are known to be exact integrals of the energy partial- 
differential equation. As has already been stated, the velocity profiles 
and stagnation-enthalpy profiles are respectively of sixth and seventh 
degree here. In reference 13, both profiles are of the more customary, 
but usually less accurate, fourth -degree type. Finally, it is of interest 
to observe that in the present analysis only a single boundary-layer 
thickness is used, instead of a dynamical and a thermal boundary-layer 
thickness, such as used in references 13 and 15 . This, however, does not 
necessarily impose any undue restrictions on the thermal profiles, since 
the latter have here been permitted to contain an additional coefficient, 
not determined in advance by the boundary conditions. This additional 
coefficient, which must be determined by the differential equations, 
replaces the thermal boundary-layer thickness as a parameter. This 
explains why the stagnation-enthalpy profiles are here of one degree 
higher than the velocity profiles. For Erandtl numbers near unity, as 
is the case for air, it has been found (cf., for example, refs. 13 , 15 , 
and 16 ) that the ratio of dynamical to thermal boundary-layer thickness 
is close to unity. Consequently, the present analysis should yield 
sufficiently accurate results for such cases . The advantages of using 
a single boundary-layer thickness are that the equations thus obtained 
are somewhat simpler and that it appears somewhat easier to solve for 
the additional coefficient than for the second boundary -layer thickness, 
especially for flow without a pressure gradient. 

This work, carried out at the Polytechnic Institute of Brooklyn 
Aeronautical Laboratories, was sponsored by and conducted with the 
financial assistance of the National Advisory Committee for Aeronautics . 

The authors are happy to acknowledge the aid of Mr. Richard Grape in 
carrying out the numerical calculations and in checking the equations. 


■4his is due primarily to the fact that in the two-thickness method, 
the ratio of the thicknesses appears in the equations in a complicated 
manner. 
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constants in equations ( 50 ) and (51) 
constants defined by equations (63) 
coefficient of "P 5 in velocity profile (eq. (20)) 
constant "average" value of &2 

coefficient of tJ in stagnation-enthalpy profile (eq. (21)) 
constant "average" value of b-j_ 

proportionality factor in temperature -viscosity relation 
(eqs. (9) and (10) ) 

constant "average" value of C 


34 


d| 


average skin-friction coefficient, i 


f x 

J (|j.du/dy) 0 dx 


( 1 / 2 ) 

local skin-friction coefficient, (p ^ u /^y) o ^(p co u 00 ^)2 


coefficient of in given wall-temperature distribution 

(eq. (Al)) 

specific heats at constant pressure and at constant volume, 
respectively 

coefficient of in equation (A 5 ) 
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F 1 ,F 2 ,F 3 integrals appearing in equations (ll) and (12) and defined 
thereafter (ef. also 6qs. (6l)) 

F^,F 2 constant "average" values of F-j_ and Fg 

f^ functions of t appearing in equation (A 3 ) and defined there 

= H 0 / H l(£)j for %r = °1 = T o/ T e (cf. also eq.* (^6)) 

Gg, C & ) ^ ( S ) # C S ) parameters defined in equations (33), (3&b), 
(^5h) , and (58) 


H 


EL 


h 

k 

L 

M 




N, 


Nu 


R. 


Pr 


q. 

R 


s 

T 


stagnation enthalpy, (u 2 / 2^) + c p T 

stagnation enthalpy corresponding to equilibrium wall 
temperature 

local heat-transfer coefficient (eq. ( 40 ) ) 
coefficient of heat conductivity 
characteristic length 
Mach number 

constants defined by equations (67a) 

Nusselt number 
Prandtl number, p.Cp/k 
static pressure 

heat-transfer rate, -(k dT/dy) 0 
gas constant, Cp - Cy 

Reynolds number based on L, p u L/ll 

Reynolds number based on x, p u x/u. 

' 00 00 / 00 

Sutherland constant; 216 ° R for air (cf . eq. (10)) 
absolute temperature- 
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t 

u,v 

**y 

7 

8,8 t 

A = R L (5 t /L) 

P 

v 

t 

P 


equilibrium wall temperature for zero heat transfer 

transformation variable, defined by equation (13) or (l4) 

velocity components in x- and y-directions, respectively 

coordinates parallel and normal to surface, respectively 

constant parameters, depending on II, and defined after 

Pr 

equation (A4) 

ratio of specific heats, Cp^c v ; 7 = 1.4 for air 

boundary-layer thicknesses in xy- and xt-planes, respec- 
tively (cf. also eqs. (15) and (l6)) 

recovery factor, defined by equation (42) 

2 


coefficient of viscosity 
kinematic viscosity, p/p 
dimensionless distance along wall, x/L 
mass density 

dimensionless variable, t/8^ 


Subscripts : 

0 values at wall; for example, T 0 

1 local ^values at outer edge of boundary layer; for example, 

and 

values in undisturbed free stream; for example, M (Also, 
cf. footnote 2) 


00 
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GENERAL EQUATIONS 


Hie following equations describe the steady, two-dimensional, laminar- 
boundary -layer flow of a compressible gas along a slightly curved wall: 

Momentum, equation in the x-direction: 


puujj. + pvuy = -p x + (jiUy) 


( 1 ) 


Momentum equation in the y-direction: 


= 0 


( 2 ) 


Continuity equation: 


(pu) x + (pv) y = 0 


(3) 


Equation of state: 


p = pRT 


00 


Energy equation: 


puc p T x + pvepTy = up x + + puy 2 (5) 

Equation (2) implies that the static pressure within the boundary layer 
is the same as the static pressure in the potential flow just outside the 
boundary layer. Thus, for slightly curved walls 


-P x = -dp/cbc = P 1 u 1 (*i 1 /dx) 


( 6 ) 


where the subscript 1 denotes local conditions at the outer edge of the 
bo unda ry layer. It may be noted that the assumption that the pressure 
in the boundary layer is equal to the potential flow pressure must be 
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modified at high Mach numbers and very low Reynolds numbers (cf. ref. 17 ). 
The methods to be used in the analysis of such flows are as yet not clear. 
As is usual in aeronautical problems the coefficients of viscosity and 
heat conductivity are considered to be known functions of the temper- 
ature, while the coefficient of specific heat at constant pressure is 
assumed constant. 


It is convenient in the method of analysis to be presented here to 
rewrite equation (5) . If equation (l) is multiplied by u and added to 
equation (5), it is readily possible to put the resulting equation in 
the form 


+ pvBy) 



(7) 


provided the Prandtl number is assumed constant. For most gases the 
actual change of Prandtl number with temperature is sufficiently Bmall 
to justify this assumption, especially since the experimental determi- 
nation of the heat conductivity appearing In the Prandtl number Is subject 
to considerable error (cf. ref. 18 ) . 

Before discussing a method of treatment of equations (l), (3), (4), 
(6), and (7) the mass density and coefficient of viscosity will be written 
explicitly in terms of the temperature. From equations (2) and (4) one 
obtains the equation 



(8) 


For mathematical convenience the viscosity-temperature relation suggested 
in references 5 and 19 is extended here to the case of pressure gradient. 
Thus, within the boundary layer 


where 



c = (V t 4 1/2 ( t « 



(9) 


(10) 


With this choice of C the approximation of equation (9) to the actual 
viscosity -temperature relation is good in the neighborhood of the wall. 
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where the skin friction and heat-transfer properties are determined, and 
is, for most problems, considered satisfactory thro u ghout the boundary 
layer. It may be pointed out that this viscosity- temperature relation 
is equivalent to using where C is similarly defined 

in terms of T 0 and Tq, if it is assumed that the ratio is 

calculated by means of the exact Sutherland relation. 

By integrating equations (l) and (7) with respect to t over the 
boundary-layer thickness t = 0 to t = and using the boundary con- 
ditions u=v=0 at t=0, together with smooth transition of the 
velocity and temperature profiles to their local main-stream values, the 
following differential-integral equations are obtained: 




r , , -il 


V + F 1 ( los e Pi)' + (lPS e u l)' 

F i + 1 1 + z-k 1 "i 2 ) f 2 


c ( p “/ Pi ) (w ui ) ( Ti / T “) (“N r 


(n) 





(12) 


where 
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A ' ( 6 t/ L ) 2 ( p „ u «, L / l ‘») 
6 = x/L 

T = V 5 t 


The prime denotes differentiation with respect to £, and t is the 
Dorodnitzyn variable defined by the transformation 


t = 




or by the inverse 


Moreover, 


or 



( 15 ) 


DA) 


( 15 ) 


(16) 


The quantities p^/p ^ l 11 equations (ll) and (12) are 

related to uq^u ro , vhich is a function of | prescribed by the potential 

flow about the body in question. Thus, in accordance with the usual 
isentropic flow relations, ^ 


O 

c In applying these relations and interpreting infinity conditions, 
the presence, in the supersonic case, of attached or detached bow waves 
should be considered. "Infinity" conditions should then be taken as con- 
ditions just behind such shocks on the streamline which in potential flow 
would form the body under consideration. 
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“l - (»i / u 4 M « 


■<i + 


(r - « ( M „ 2 / 2 ) 


1 - (V 2 -/ 


^1 - 1 / 2 


( 18 ) 


Approximate solutions to equations (ll) and (12 ) can be found by- 
assuming u/u-^ and H/H-, as functions of t which satisfy certain 

boundary- conditions at r = 0 and r = 1 to be discussed in detail 
later. In these profiles one parameter in addition to Tv is left unde- 
termined after the boundary conditions are satisfied, and these two param- 
eters are considered as the two unknown functions of £ to be determined, 
by equations (ll) and (12) . 

In setting up the final differential equations, examination of equa- 
tions (ll) and (12) indicates that, in addition to selecting the profiles 
for u/Up and H/H^ with their appropriate boundary conditions, one 

must prescribe (|), M m , and Np^,. Moreover, the values of the 

fluid properties at infinity must also be specified. From the solutions 
to equations (ll) and (12) all quantities in the £ T-plane can easily be 
found; the transformation, equation (l4), can be applied to obtain the 
corresponding quantities in the xy or physical plane. For this purpose 
equation (l4) may be conveniently written as 


y/L = (8 t /L) < (H/Hj 1 + (7 - 1) {vi^j 2 ) - (7 - 1) (M^/ 2 ) (u/u^ 2 l dr 

l J(l 9 ) 


The boundary conditions that the velocity and stagnation-enthalpy 
profiles are to satisfy will now be discussed. As indicated in the 
"Introduction" a sixth-degree polynomial in t for the velocity and a 
seventh-degree polynomial for the stagnation enthalpy are suggested on 
the basis of the critical study of the use of integral methods in 
compressible-boundary-layer analysis presented in reference 1. Thus, 
let 
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U / U l = 51 a t 3 T ' 5 

' J=o 


( 20 ) 



( 21 ) 


and b coefficients except one are determined from 
J 


All of the a. 

d 

boundary conditions at r = 0 and t = 1. "In setting up the differential- 
integral equations (ll) and (12 ) , the boundary conditions -which must be 
specified to obtain a unique and exact solution to the governing partial- 
differential equations have been applied. In the £ t - plane these are: 

At t = 0, 


u/ui = vj\i = 0 
S/Si - (H 0 /B^)(|) 


'■'N 




(22) 




and at t = 1, 


u / u l = H /^l = 1 

( u /S> T ■ - 0 








(23) 


The boundary conditions specified here are the usual ones in heat-transfer 
problems. It is possible to specify instead of either J 

or some relation between these two quantities as in the case of radiation 
cooling. In principle there is no loss in generality here, since 

is left arbitrary and nothing is specified with respect to 


mi 
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In addition to requiring that the assumed profiles satisfy boundary 
conditions (22) and (23), it is customary in the integral methods (cf., 
e.g., ref. 1) to select the profiles so that at the wall and at the edge 
of the bo undar y layer seme of their derivatives have the same values as 
those which an exact solution to the orig in al partial-differential equa- 
tions under boundary conditions (22) and (23) would yield. For the pro- 
files of equations (20) and (21) these additional conditions can be 
obtained from equations (l) and (7) and from differentiation of these 
equations with respect to t. In this way one obtains the following 
boundary conditions to supplement those of equations (22) and (23): 

At r = 0, 


c (V T 4 - - A ( p i/ p 4 (ys) £ + (r - 1)M i 2 /?] ( u i/ u ») ' (2k} 


(V%) W m - w„w T 


(25) 


1 + (7 - 1) 





("IN* H FT ( p l/ p „) ( u A) t PoA) ’ 



3(i - y (r - i)m l 2 (u/x l ) t (u/ U i ) 


1 + (7 



(27) 


At T = 1, 



PA) 


= 0 


TTT 


(28) 
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With, equations (22) to (28) f ourteen of the fifteen coefficients in 
the velocity and stagnation-enthalpy profiles of equations (20) and (21) 
can be determined in terms of A, of the fifteenth coefficient (to be 
chosen as b-j_), of the prescribed functions ( u i/u<»)(0 sonh (6 ) > 

and of the prescribed parameters M^, and Npp. 

As explained in the "Introduction, " it will be observed that only 
one botindary-layer thickness is here assumed and that thus no 

distinction is made between a velocity and a temperature boundary-layer 
thickness. This, however, does not impose an undue restriction on the 
shape of the temperature profiles, since the latter (see eq. (21)) are 
here permitted to retain an additional parameter to be determined, not 
by any boundary conditions, but by the differential equations (viz, 
eqs. (ll) and (12)). 

Examination of equations (24) to (28) indicates that considerable 
simplification is achieved in the boundary conditions, and therefore in 
the analysis, under two special cases: (a) The flow without axial pres- 

sure gradient ^u^/u^ = lj , which also corresponds to the practically 

interesting case of the supersonic flow over a thin wedge, as well as 
to the subsonic (or supersonic) flow over a flat plate and (b) the air- 
foil with pressure gradient, and constant wall temperature, in a fluid 
with Prandtl number equal to unity. Although the assumption of constant 
wall temperature does not appear to be exactly realized on practical 
wings, it appears desirable to investigate the boundary-layer character- 
istics under these assumptions since the pressure-gradient, heat-transfer 
case in a compressible flow has not yet been analyzed entirely satisfacto- 
rily from the standpoint of simultaneous accuracy and simplicity. 


FLOW WITHOUT AH AXIAL PRESSURE GRADIENT 


In this section a general solution based on the previously derived 
integral-differential equations will be obtained for the special case of 
flow without axial pressure gradient. The accuracy of this solution is 
then investigated by comparing it with more exact solutions. 


General Case 

In the case of the flow without axial pressure gradient and with 
general but constant Prandtl number and variable wall temperature the' 
integral-differential equations (11) and (12) become: 
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(*]»*’ + V X - C ( U / U l) 


( F 3^) V + * 3 '* - ( C / N Rr) ( E /=l) T 


Here M x s Moo, T-, = T^, and so forth. The subscript » will be used 
to denote the free-stream conditions for this case. 

It is convenient to consider A and b^ as the unknowns to be 

determined by equations ( 29 ) and (30) . Then the velocity profile satis- 
fying boundary conditions (22) to ( 28 ) is 

u/u x = 2r - 5 t^ + 6t^ - 2t^ (31) 

while the bj coefficients of the stagnation-enthalpy profile become: 
b o = Hq/Hco = G 1 


b 2 = 2G 2 


bj = 2G^A 

b^ = 35(l - Gq) - 20G 2 - 8 GjA - 20b x 
b^ = -84 ^1 — G^ + 40G 2 + 1 2GjA + h-^b^ 
bg = 70^1 - - 30G 2 - 8 G 5 A - 36 b x 

b^ = —20 ^1 — Gqj[ + 8 G 2 + 2G^A + 10b x 



1 6 


NACA TN 3157 


where the G's are known functions of | defined by: 


°1 " ( H q/ H ») <*> 

g 2 = (! - H Pr) (7 - i)^ 2 / 1 + (7 - l)(M M 2 /2)j 
G3 = (1/6)1?^ G 1 '/c 




(53) 


From the definition of F-^ and from equation (31) it is found that 


F x = 985/9,009 


(34) 


Similarly, 


F 3 = (31/126)^1 - G^ - (953/l80,l80)G 5 A - 

(302/9,009)G 2 - (821/12, 0'^b-L (35) 

With the value of F-]_ given by equation (34), A(|) can be found 

directly from equation (29). Thus with the requirement that A(0) = 0, 

X = 40^ (36) 


where^ 


c ± = c de 
Jo 

Equation (30), which gives b 1 (i), becomes: 

f 3 1 - (Y™Rr) (“/OlK + 5/ a )(7 c i) ■ 0 (5T) 

5lt may be noted that in this analysis the temperature-viscosity 
proportionality factor C, unlike in reference 5, can be permitted to 
vary along the plate (i.e., C = C(|)) without any mathematical difficulty. 
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The solution of equation (37) satisfying the requirement that b-j_ he 
finite at 5 r 0 is 


-Pn pi Pi 

b l = C 1 / C 1 
J 0 


(38a) 


where 


P x = ( 1 / 2 ) + ( 985 / 2 , 463 ) (i/n^) 

0^(6) = (12,012/821) ^(31/252) (l - Gi) (G/C-jJ - 

G 1 '|^(3l/l26) + (953N pr /29,550) (3/2) - (c^ C 2 j j> - 
(953N pr /29,55°) (Ci/c)G 1 n - (151/9,009) (c/cj 




(38b) 




From. equation ( 38 a) with any prescribed temperature distribution 
G-^f;) at the wall, b-j^ = b-^l) can be readily determined. From equa- 
tions (32) the remaining bj coefficients can then be found as functions 
of | and of the prescribed flow parameters, for example, N-^ and M^. 
With the aj and bj coefficients and A all determined, the velocity 

and stagnation-enthalpy distributions throughout the |T-plane are known. 
By application of equation (19), these distributions in the xy or phys- 
ical plane are easily determined. 

5he average skin-friction coefficient Cj for the length L can be 
calculated from equations (31) and (3 6 ) as 


C f 



1.322 




(39) 


where 



C d|. 
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The local heat-transfer coefficient h is usually defined as 

h= -(kTy) y(T 0 - T e ) (40) 

Using solution (3 6 ) for A, equation (40) yields: 



; where q is the recovery factor defined by the equation 

I 

I 

T e/ T l = 1 + ‘ = ( H e/ S l) f + (7 “ 1 )( M 1 2 / 2 ) (^ 2 ) 

The (nondimensional) Nusselt number is thus 



The value of T) .follows by setting b^ = 0 and solving for the 
constant wall temperature, as determined by G]_* Equations (33 ) , (38), 
and (42) with 7 = 1.4 are thus found to imply: 


71 = 1 - 0.272(1 - 


( 44 ) 
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It should be pointed out that equation (44) yields for Np^ = O .72 

a value of t\ = 0.924, -which must be compared with the value of 0.845 
given by the mathematically more exact method of Chapman and Rubesin. 
Although this result must be considered a consequence of the appro xima te 
method of analysis used here, calculations involving rates of heat trans- 
fer (such as Nusselt number) for given values of the wall temperature 
will not be directly affected by the equilibrium wall temperature and, 
therefore, will not be appreciably affected by any ina ccuracy in tj.4 
This is borne out by the numerical example and the comparison with an 

exact solution, given in a subsequent section here, with = 0.72. 

J Pr 1 

It may nevertheless be convenient for certain applications to express 
the wall temperature in terms of the predicted equilibrium temperature. 
With the additional simplification of replacing variable C by an average 
value C equation (43) becomes: 


where 


%u/^ - 0 - 2 97 ^[l - (T c /T e ) \ Pl jT 1 05 ( 45 a) 


°5 = I 


Pq-1 


|l - (T 0 /T e )J “ [2 + ( 60 , 039/152, 675)Np^ (T 0 / /T e )'5 - 


"e 2 


(4o,026/i52,675)Rpr(T 0 /T e ) "S 


(45b) 


In deriving equations (45a) and (45b), the following relation between 
and T 0 /T e (cf. eq. (42)) was used: 


°1 - ( T o/ T e) 


1 + 



(7 - 1) (Mp 2 ^) 


(46) 


^Because of the inaccuracy in q, it will be found, of course, that 
when the wall temperature has the same value as the true equilibri um 
temperature, the present equations will not lead exactly to a zero heat 
transfer at the wall. However, in such a case, it will be found that the 
calculated rate of heat transfer will be relatively small . 
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It should be observed that the recovery factor T) does not appear in 
solution (^5a), although T Q/ /T e is prescribed. 

If it is desired to calculate the actual rate of heat transfer q 
(in British thermal units per second per square foot) for a given 
ratio T Q /T e from equation (45a), then the Nusselt number as determined 
by equation (45a) should be multiplied by - T e j , where T e Is 

the actual equilibrium wall temperature (as given, e . g. , by exact theory 
or by experiment). In practice, however, the wall temperature T Q (or 

the ratio T 0 /T ro ), rather than the ratio of wall temperature to equilib- 
rium temperature, will usually be prescribed. In that case, the rate of 
heat transfer can be directly calculated (without any intermediate use 
at all of the recovery factor q) from the expression: 


* a -(* *W) 0 = -C (1 ♦ ifi | ^ (47) 


The quantity b-^ can be calculated from equations ( 38 a) and ( 38 b) for 
flow over a flat plate, with the observation that 


°i - (Vi) i 1 + H 1 “i 2 ) 1 (48) 


Equation (47) has been derived with the use of the following relation 
between temperature and stagnation-enthalpy profiles, resulting from the 
definition of the stagnation enthalpy H: 



In order to determine the velocity and temperature profiles in the 
physical (xy) plane, the inverse transformation given by equation ( 19 ) 
must be used to find the coordinate y associated with the variables 5 
and t, while the x-coordinate may be found directly from the definition 
of |. 
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If the wall -temperature distribution is prescribed as a polynomial 
In the distance | along the wall, then the Nusselt number and temper- 
ature profiles, following from equations (45a), (45b), and (49), can be 
directly calculated, without any quadratures, from the explicit expres- 
sions given in the appendix. 

It is of interest to observe5 the wall -temperature distribution 
required for zero heat transfer according to equations (45a) and (45b). 
Thus, by setting = 0, a second-order 1 Inear differential equation 

in T Q /T e is obtained, whose general solution is: 


- 2 . = 1 + A£ 
T e 


- 1/2 


+ B| 


-7-63/Nj* 


(50) 


where A and B are arbitrary constants. Equation (50) implies that 
zero heat-transfer conditions will be satisfied not only by a uniform 
wall temperature T e but also by a wall temperature varying inversely 

as a certain power of the distance £ along the flow. In particular, 

the heat-transfer rate vanishes for 1 - (t o /T 6 ) « £ (regardless of 

the Prandtl number). This result was also obtained by Levy (ref. 20) for 
incompressible flow, using an analysis restricted to power-law wall- 
temperature distributions. The wall temperature according to this distri- 
bution would be infinite at £ = 0, but Eckert (ref. 21) has reasoned 
that a heat-line source or sink placed at the region where | = 0 might 
result in a temperature distribution of this type. The case of upstream 
cooling such as analyzed in reference 22 Is the physical counterpart of 
the abstraction of this line source or sink. Equation ( 50 ) indicates 
that there is another type of wall -temperature distribution for zero heat 

-7.63/Npj. 

transfer, namely, 1 - (T Q /T e J « £ . This is not included in 

reference 20, where negative powers of £ for T 0 (£) lower in absolute 
value than 7*65/%^, (approximately 10 for air) were assumed. It should 

be noted that the results obtained in the present analysis are valid for 
compressible flow (i.e., for any Mach number). Finally, it may be noted 
that temperature distributions of the form of equation ( 50 ) are not 
included In polynomial distributions, such as those to which reference 5 
is restricted. 


^The observations in this paragraph were pointed out to the authors 
by Dr. M. Blocm. 
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Constant Wall Temperature and Np^. = 1 

A special simple case of the above solution for flow over a flat 
plate is of interest here. This is the case of constant wall temperature 
and Prandtl number equal to unity. This case will now be shown to give 
a particularly simple check on the accuracy of the method presented here. 

It is well known that if the Prandtl number is assumed equal to unity 
then equations (l) and ( 7 ) imply for the flat-plate case that 


H = A + Bu (51) 

where A and B are arbitrary constants. Equation .(51) is an exact 
stagnation-enthalpy relation which satisfies, the boundary conditions of 
constant wall temperature both with and without heat transfer. Therefore, 
one test of the approximate method presented^ above Is its automatic 
reduction to this exact relation, if the conditions under which this 
relation is known to be exact are imposed. This test will now be applied. 

For the case of constant wall tenperature and a Prandtl number of 
unity, equations ( 38 b) and ( 33 ) imply: 

04 = (4,433/2,463) (1 - aJd/O 

Hence, equation ( 38 a) yields: 

b-L = 2(l - G 1 ) (52) 

Furthermore, from equations (32), 

bg = bj = 0 

b 4 - -5(1 - 03.) 

bj = 6(1 - GjJ 
b 6 . -2(1 - O3) 
by = 0 
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When these b j coefficients are substituted into the stagnation-enthalpy 
profile, one obtains 

H/% = G ± + (l - ( 2 t - 5 t 4 + 6x5 - 2 t 6 ) ( 5 A) 

and thus, by comparison with equation ( 31 ) , 

H/H; 1 = G L + (l - G x ) (uju-^ (55) 

Equation (55) is exactly equation (51) in. nondimens ional form, with 
the boundary conditions that H = H Q when u = 0 and H = H-^ when 

u = Uj_. It is thus seen that the reduction to the exact energy equation 

for this special case is automatic and rigorous. It has already been 
shown in reference 1 that the agreement between the results of the inte- 
gral method using a sixth-degree velocity profile in conjunction with 
equation ( 55 ) (i.e., for constant wall temperature and Prandtl number 
equal to unity) and the results of the exact solution of reference 5 is 
good. The skin friction and heat transfer are for all practical purposes 
exactly predicted and the laminar -boundary- layer stability limits are in 
good agreement. It can therefore be expected that the more general inte- 
gral method presented here, valid for cases of heat transfer in a pres- 
sure gradient, will be accurate in more general cases, since the two 
integral methods are in complete agreement for the special case just 
considered. 


Comparison With the More Exact Method of Chapman and Rubesin 
(Variable Wall Temperature and Npp j 1) 

In reference 5> Chapman and Rubesin have applied their analysis to 
the case of a variable, prescribed wall temperature. To simplify the 
numerical analysis the proportionality constant C was assumed by them 
to be equal to unity. The Rrandtl number was assumed to be O. 72 . In 
the notation used here the wall temperature was expressed there as 


T c /T e = 1.25 - 0.831 + 0.331 2 


(56) 


Thus for equation (Al) (see appendix) c Q = 1.25, c-^ = -O. 83 , and 

c -2 = 0.35. The Russelt number variation with | given by equation (A2) 

is shown in figure 1, along with the results of the more exact Chapman- 
Rubesin calculation. It can be seen that excellent agreement is realized. 
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It will be noted from equation (39) with = 1 that the average skin 
friction obtained by the present equations is in excellent agreement 


with the exact value 



= 1 . 328 , given in reference 5* 


In comparing the velocity and temperature profiles predicted by the 
two methods it must be kept in mind that for strict comparison the wall- 
temperature distribution should be precisely the same in the two methods. 
However, since Chapman and Rubesin assumed only the ratio of wall temper- 
ature to equilibrium temperature and since the two methods do not give 
precisely the same recovery factor, and hence not the same equilibrium 
temperature, it is necessary for this comparison to express the wall 
temperature assumed by Chapman and Rubesin, as well as that used in the 
present report, in -terms of the free-stream temperature. 

In figures 2 and 3 the comparisons of the temperature and velocity 
profiles predicted by the two methods are shown for three values of | 
and for two Mach numbers . It can be seen that excellent agreement is 
realized. 


It may therefore be concluded from these numerical results that the 
momentum-and-energy integral method given here will yield good results 
with little computational difficulty when applied to the case of zero 
axial pressure gradient, variable wall temperature, and general but 
constant Rrandtl number of the order of magnitude of unity. 

From a practical point of view, it should be kept in mind that the 
solutions developed here are based on the viscosity- temperature rela- 
tions ( 9 ) and (10), which are an approximation to the actual relation 
for air. Because of relations ( 9 ) and (10), the results obtained here, 
namely equations (39) and (^5a), indicate that, for a fixed wall temper- 
ature, the skin-friction coefficient and the Nusselt number will be 
independent of Mach number. For the Sutherland viscosity-temperature 
relation, however, this will not be quite valid (cf. ref. 8). 


FLOW WITH PRESSURE GRADIENT 


For flow over a curved surface, such as an airfoil, where the local 
velocity distribution u^y/u^ outside the boundary layer is not constant, 

but may be considered as a given function of the distance x along the 
wall, the ordinary differential equations (ll) and (12), unl ike the case 
of zero pressure gradient, can no longer be easily solved for general 
Rrandtl number and arbitrary distribution of temperature along the wall. 
For such flows, however, considerable mathematical simplifications occur 
when the Prandtl number is unity and the wall temperature is uniform 
(cf. especially eqs. ( 26 ) and ( 27 )). 
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The assumption Rpp = 1 has the additional theoretical advantage 
that, according to equation (Ml-), it automatically leads to the known 
exact value of the recovery factor, namely, tj = 1, for this value of 
the Prandtl number. In fact, the parameter G^ now ass ume s particular 

physical significance, since for Npj. = 1 equation ( 46 ) implies: 

Gl - T 0 / T b (57) 

Consequently, the case of Npj. = 1 and uniform wall temperature will he 

treated now in detail. It will he shown that for such a case a relatively 
simple approximate solution of the equations can be obtained for any 
given velocity distribution outside the boundary layer. 


General Approximate Solution for 
Npr = = Constant 

For this case, with (l) arbitrary, the coefficients a^ 

and bj in equations (20) and (21) for the velocity and stagnation- 
enthalpy profiles can, by virtue of boundary conditions (22) to (28), 
all be expressed in terms of a 2 and bq, where b 1 remains arbitrary, 

while a£ is given by 


2-7 

a 2 = -C 1 /^) Hq^A = G 6 (g)A (58) 


where Gg(| ) is a given fi unction'. The profiles can then be given in 
terms of and bq as follows: 

u/uq = (zr - + 6t^ - 2t ^ _2t + 5^ - IOt^ + IOt^ - + 


jb^fi Gq^-T + 10 t 5 - 20t^ + 15 t 5 - 4t 6 ) 


(59) 
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=/®L ' »i + (1 - Oj) (?5 t 4 


b^ t - 20t m ’ + 45 t^ 


84t 5 + 70 t 6 - 20r^j + 
36t 6 + 10t7) 


(60) 


With expressions (59) and (60), the following explicit expressions 
for F^_, f 2 , and are obtained: 


F, = 


F o = 


F_ 


0.1093 + 0.00211a 2 - 0.000622a 2 2 +• 0.000412^-^2^ - 
0.0000095^^2^ 2 - 0. 000155 (t>2. a 2 2 / G l) 

0.395 - 0.500^1 - + 0.107b 1 + 0 . 0212a 2 - 0.00062a £ + 

0.0028^^^ - 0.00015 - 0.0000095^^ 

- C 1 - G i)[ c 


> ( 61 ) 


0.246 - 0.015a 2 - O.OOlSl^ag^ - 


0.0683 - 0.00324a 2 - 0.00041^^ 


With expressions (59) ^ (60) , and (6l) inserted in equations (ll) 
and (12), two ordinary differential equations for X(|) and b^(£) axe 


obtained. Although these can be solved numerically for a given 


W 6 


the process may be tedious . A relatively simple general approximate 
solution of these equations will, therefore, be derived. 


Equation (ll) can be solved approximately for A by assuming that 
Fq_ and F 2 can be replaced there by constant "average" values Tj^ 

and Fg over the distance | . This is justified by the fact that the 
variable terms there, which are proportional to a 2 and b^, are rela- 
tively small (cf. eqs. ( 6 1)). This is equivalent to replacing a g 
and b-j_ by constant "average" values and b-j_ for this purpose. 
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With, equation ( 59 ) for the velocity profile and equation (58) for a^, 

equation (ll) then becomes the following linear ordinary dif ferential 
equation in A: 


( F i/ 2 ) v + m f i( p i7 p i) + (v/4 


”1 + H 1 - 4 - 


2C 


(Pco/P 1) ( T l/ T oo) 4/4 


(62) 




(63) 


where cp^ is a constant, given by: 

”1 - 5 1 + ? 2 - 5°1 - N 60 ) 

= A + Db x 

A = 0.3G 1 + 0.00*1-38 + 0.0232a 2 - 0.0012*l-a^ 2 
D = 0.0905 + (0.0838 - 0 . 004 58^ (& 2 j30G^j 


With relations (17) and (l 8 ), the solution of equation (62) satisfying 
the condition A = 0 or finite (if u-j_ = 0 at £ = 0 ) at the leading 

edge (| = 0 ) is found to be: 


A = 


, h “-i; 

tej, k w m 

\i y ^ ni.h 

F-, 1 7-1 V 

(44 (44 1 


d£ 


m 
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Equation (64) is similar in form to equations obtained for zero heat 
transfer in references 3 and 23 and for heat transfer, but with fourth- 
degree profiles and two boundary-layer thicknesses, in references 13 
and 15 . 

To obtain a general approximate solution for b^(|), V can first 

be eliminated from equations (ll) and (12) . Assuming F^' = F ^ 1 =0 

for this purpose and using equation ( 58 } , the following equation is thus 
obtained : 



Noting, in advance, that b^ will usually have a value approximately 
equal to 2^1 - G^, “b^ can be replaced by 2^1 - Gq^b^ in the 
relatively small term 0 . 0000095 ^i a 2 /^l) ^ appearing in F-^ and Fq 
and in the small term 0.0004 lb^ag/G^^ appearing in F^. Expressions (6l) 
for F-^, Fg, and F^ can then be written as linear functions of b-^, 
with coefficients as functions of Gp and a^. Ey substituting these 
expressions into equation ( 65 ), the following quadratic equation in b^ 
is obtained: 

mb-j^ + jb^ -1=0 (66) 

where 


2a 2 1 a -1 *' 

m ■= x 10 


65.90 - 4.496a P + 2.705 Si . 0.2096 + O.Ol452/^\ + O.OO56 ^ - - 

G i. G i . ,, vhy % 

0.00052 ^ - 0. ooois 

G X 2 \flj j 

3 = 0.24602 + 0.07917ag - 0.005868a 2 2 + 10 ^284.5 + 42.26a 2 - 2.067a 2 2 + 

0.0465a 2 3 - 7.906 — + 0.9086 Si- + O.ChS^/^] - 0.0232 _ 0.00284 

°i °i \<h) . Qi 

(1 - 0.24602 - 0.0l496agj 1 + ^-0.10495 + O^ + 0.021l6a 2 - 0.00062l6a 2 2 j 


> (67a) 


l =• 2 
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If terms which will ordinarily be negligible are rejected, equations (67a) 
can be slightly reduced, thus: 


2ap 

m = X 10 

G-, 


63.90 - 4 . 496 a 2 + (2.705 - 0.2096a 2 ) 


, . 2ap _4 

J = 0.24602 + a 2 ( 0.07917 - 0.005868a^ + — X 10 X 


284.5 i- (42.26 - 2 . 067 a 2 )a 2 - ( 7.906 - 0 . 9086 a 2 ) ^ 
l = 2(l - Gq) (0.24 602 - 0.0l496a 2 ) ^-0.10495 + 0.3G 1 

a 2 (o.0211o - 0.00062l6a 2 )Jj> 

Hie solution of equation ( 66 ) i 6 : 


> 


(67b) 


b = - 
L 2m 


(±) 2 + 1 

i2m / m 


1/2 


( 68 ) 


(The physically appropriate root will in general be that which is closer 
to the value 2 (l - Gq) . j 

Equations ( 64 ) and (68) represent a simple approximate solution of 
equations (ll) and (12) for Np^, = 1 and uniform wall temperature 

(Gq = Constant) . Their application will involve at most numerical 
integration . 

In applying equation ( 64 ), a reasonable "average" constant value 
for a 2 , for any given (uq/u^j (| ) , Gq, and Meo, can usually be obtained 

by considering equation (58) for a^C/GqA and equation ( 64 ) for ~h/C. 

An "average" value bq for b l' to be used in evaluating cp^ (eq. (63) ), 
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can be obtained by considering equation (68) for b]_. For objects with 

sharp leading edges, for which A = 0 at |_= 0, it will ordinarily be 
found that b^ «* 2^1 - G-^ very roughly. After A(|) has been found 

from equation (64), b^(l) can be easily determined frcm equation (68).^ 

With l(|) and b-j_(|) determined, the boundary-layer characteristics 
can all be straightforwardly calculated. The local skin friction will be: 


(°/\/*)(VM N u “) R i' 1/a 

( 69 ) 

The Nusselt number, giving heat-transfer properties at the wall, will be 


(p8u/8y) 

Cf 1 h a. 2. = 4 

1 1 2 
« Poo U co 


" ( a s/5) - (V2/ 60G l) 


%u = 


(k8T/8y) Q L 
^00 (T 0 “ ^e) 



(70) 


The velocity and temperature profiles follow from equations (59) and ( 60 ), 
in conjunction with equation (49) for T/Tp, and equation ( 19 ) for trans- 
forming to the physical plane. 

In accordance with conclusions reached in references 1, 2, and 3> 
where the methods applied were essentially the same as in the present 
investigation (viz, sixth-degree profiles, together with a solution of the 
form of eq. (64)), it may be expected that the results obtained by using 
the equations developed in this section will be sufficiently accurate for 
most practical purposes, including stability calculations. However, for 
flow over a blunt-nosed object (i.e., "stagnation flow") or for determi- 
nation of the separation point in an adverse pressure gradient, these equa- 
tions may have to be modified to yield still greater accuracy (cf . ref. 3, 
where such modifications are shown for zero heat transfer). For stagna- 
tion flow (uWu^ = k|\ , in fact, fourth-degree profiles give satisfactory 


f) 

In equations ( 67 a) and ( 67 b), the exact expression ( 58 ) for 
should be used and not any "average" value. 


a 2 (0 
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accuracy (ref. 3). Consequently, stagnation flows with, heat transfer 
could he calculated by the method of reference 15 , where such profiles 
are used for this case. 

It should, finally, be observed that equations (6^) and (68) have 
been based on the assumption that Fq ? F 2 , and Fj are approximately 

constant, or vary slowly, along the flow. Regarding F-j_ and F^, this 

can be interpreted physically as an assumption that the ratio of the 
momentum thickness to the boundary-layer thickness, and the ratio of the 
thermal momentum thickness to the boundary-layer thickness, will be 
approximately constant along the flow in the xt-plane. This approximation 
has been justified by the relatively small variable a 2 terms (replaced 

here by terms in a g) in expressions (6l) . In ordinary cases, the a 2 

terms will actually be found to be sufficiently small, either individually 
or collectively, so that the approximating assumptions made here will be 
satisfied. In any doubtful case, however, one can in general easily 
check, a posteriori, the validity of the approximate solutions developed 
here by computing the variable a 2 terms with the approximate solution 

obtained and ascertaining whether these terms have a sufficiently small 
net effect, relative to the other terms, to be replaced by constant 
"average" values. Assumptions corresponding essentially to the use of 
constant "average" values of F]_ and F 2 have already been applied in 

previous analyses (refs. 3 and 15 ) • 


CONCLUSIONS 


From the analysis presented herein for the laminar compressible 
boundary layer with heat transfer in flow with and without an axial 
pressure gradient, the following conclusions can be drawn: 

1. The boundary-layer characteristics for the flow over a flat plate 
can be easily determined with sufficient accuracy from the equations 
developed here, for a given constant Rrandtl number (of the order of mag- 
nitude of unity) and a given wall-temperature distribution. 

2. For flow with a pressure gradient, the boundary-layer character- 
istics can also be easily determined with sufficient accuracy by the 
equations developed here, provided the Rrandtl number is unity and the 
wall temperature Is uniform. Here, the velocity distribution outside the 
boundary layer and the free-stream Mach number are • considered as prescribed. 
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Some modifications of the equations may, for greater accuracy, be neces 
sary for flow over a blunt-nosed object and for determination of any 
possible separation point. 


Polytechnic Institute of Brooklyn, 

Brooklyn, N. Y., December 11, 1952. 
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APPENDIX 

POLYNOMIAL WALi-EEMFERAQAERE DISHUBUTION 
If the ratio T 0 ^T e is represented as a polynomial in £, that 

is if 


T 0 /T e = 


M 


j=0 


C ^‘ 


(Al) 


then equation (4 5a) yields: 



M -L P 

^ Cjl 3 (j + Pi)' -U+ 2J + 
J— 1 L f 


(60,039/152, 675) ^[l + 13,342 (j - 1 ) 720 , 013 ] 


(A2) 


With the wall temperature specified as in equation (Al), equation (19) 
together with the use of equations (31), (32), (36), and (49) leads to 
the following explicit expression for y(|,-r): 


(y/2x) = ^9,009^/985 |jV - f x + (2f 2 /p L ] + 

(r - 1 ) (m x 2 / 2 ) h 2 f 2 - f 3 + p 3 f 4 ) + 


1 + (7 - 1 ) ( M 1 2 /^) n 



[2f2/(j + Pi)] t 1 + & 


j + 



(A3) 
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where 


f x = T 1 - + 14 t 5 - 10t 6 + ( 5/2 )t 7| 

= (4 , 433/9* 852 )t 2 [2 - l6r^ + 30 t^ - (144/7) + 5t^J 
fj = -(t 3 /i 8) [-24 + 126t2 - 192t3 + (828/7) T 1 * - 27t 5 - 50t 6 + 108tT - 

(1,008/11)t 8 + 36t 9 - (72/13)t 10 ] 

f 4 = i^/i) [2 - 12t 2 + 20t 5 - (90/7)^ + 3 t 5 J 

f 5 = (t^/ 20) [lO - 32 t + 40 t 2 - (160/7 )t 5 + 5t^J 

For convenience, the f functions and their first derivatives with respect 
to time, denoted by f ^ ' , have been calculated and are plotted in 

figures 4 and 5- 

The temperature profiles in this case can be determined from the 
explicit -expression: 


- £ - v + ^ 2 'AOj + - i >( m i 2 ^) (w - f 3 + p 3 v) 

j^l + (7 - 1)^/2)^ e^-jfi’ + - 


[2f 2 '/(3 + Pjjj jl + + (40,026/152,675) j 2 j- 


(A4) 
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where 


Pi 

02 

P 3 

P5 


(1/2) + (983/2,463) (l/Npr) 

1 - ( 1 , 208 / 4 , 433) (1 - H Er)l ( 2 Al) 

2 (l - %r) 

2^/3^ = ( 6 ^ oo6 /985)N pr 
2 + (20,013/152, 675 )Np r 


The velocity profiles may be found from equations (31) and (A3). 

It may be noted that, if desired, the ratio T 0 ^T-^ may be prescribed 
as a polynomial in |, thus: 


11 -S 

To/ T 1 = 

1 3=0 3 


(A5) 


As can be seen by comparison of equations (Al) and (A 5 ) in connection 
with equations (46) and (48), the results given by equations (A3) and (a 4) 
can still be used, provided that one substitutes 



i[ 2 


1 + (7 - 



(A6) 


there . 
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Figure 1.- Nusselt number distribution % U /^R^ versus £ (eq. ( A2 ) ) . 
Flat plate; variable wall temperature; Npp = 0.72. 
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Figure 2.- Comparison of temperature profiles in ‘boundary layer. 
Variable wall temperature; Np r = 0.72. 





Figure 3»- Comparison of velocity profiles in boundary layer. Variable 

wall temperature; Npj. = 0 . 72 . 
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